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Abstract 

We explore the possibility of using a double-tip STM to probe the single electron 
Green function of a sample surface, and describe a few important applications: (1) 
Probing constant energy surfaces in k-space by ballistic transport; (2) Measuring scat- 
tering phase shifts of defects; (3) Observing the transition from ballistic to diffusive 
transport to localization; and (4) Measuring inelastic mean free paths. 
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The single-tip scanning tunneling microscopy (STM) has become a major tool for surface 
analysis [|TJ. However, its use has been limited to probing static properties of electronic 
systems such as the local density of states on or near sample surfaces. In addition, it lacks the 
k resolution to enable one to determine the energy dispersion in band structures. Transport 
properties are also out of the reach of the single-tip STM, except for the BEEM configuration 
used to probe ballistic transport across a metal film ||. In this letter, we explore the 
possibility of realizing a double-tip STM, and describe a few important applications. 

In a typical double-tip experiment, electrons are emitted from one tip, and propagate 
through the sample, some of which are picked up by the other tip. Naturally, the propaga- 
tor or the Green function of an electron in the sample is involved. Since all single particle 
properties of the system can be derived from the Green function, one expects that a lot 
more information about the sample surface and nearby region can be learned from an ex- 
periment using a pair of tips than from using a single tip. Some important applications are: 
(1) deducing useful information about the band structure of surface states, (2) measuring 
scattering phase shifts of surface defects, (3) Observing transition from ballistic to diffusive 
transport to localization, and (4) Measuring inelastic mean free paths. 

Shown in Fig. 1 is a schematic experimental set up. The sample is assumed to be large 
enough to have a well-defined chemical potential /xq- Voltages V\ and Vi are applied to the 
tips relative to the sample, and electric currents 1% and I2 from the tips to the sample are 
measured. Like the ballistic energy electron microscope (BEEM), this is a three terminal set 
up 0. Unlike the BEEM, here the tip-sample separations, the tip biases (VI and V2), and 
the tip locations (ri and r 2 ) are controlled independently. Direct junction conductances at 
i - ! and r 2 are defined as <jj = J^l, and are given to leading order in tunneling rate as 
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where p(rj,//j) is the local density of states of the sample at the chemical potential of tip-i, 
and Tj describes tip-sample couplings as well as the density of states of the tips (which 
are routinely measured in single-tip STM experiments). Up to order Til^, there are also 
processes of coherent tunneling of an electron from one tip to the other through the sample, 
which can be measured through the transconductance defined by cr 21 = J^. The transition 
rates can easily be accounted for by the Fermi golden rule using second order transition 
matrix elements ||, yielding 

2ire 2 

^2i = r 1 r 2 ^|G(r 1 ,r 2 ;6 = yUl )| 2 , (2) 

where G(ri, r 2 ; e) is the retarded Green function of the sample, for non-interacting electrons 
at zero temperature. 

In the presence of electron-electron and/or electron-phonon interactions in the sample, the 
Fermi golden rule gives the same result if fii = /x 2 and if the sample is nonsuperconducting. 
For fii ^ /i 2 , four-point Green functions of the sample are involved to account for the 
inelastic scattering of the tunneling electrons ||. For a superconducting sample, Andreev 
processes can also contribute to the transconductance within the same order |J. With the 
direct conductance and transconductance measured at \i\ = \i 2 = f-i, one can thus obtain 
information about the retarded Green function as shown below 

\G(r 1 ,r 2 ;e = /i)| 2 = 2 ™ ° 2X p(n,//)p(r 2 ,ft). (3) 

n o x o 2 

Since constant-current-STM-images trace constant contours of local density of states, the 
last two factors can be treated as a normalization constant. 

We next identify experimental parameters from which such a measurement can be real- 
ized. As a second order process, it is clear that a 2 \ is a weak signal to detect. Taking BEEM 
as a reference, for a total tunneling current of 1-10 nA, the detecting limit is about 0.1-1 pA 
for ballistically transported electrons, corresponding to a factor of 10~ 3 — 10~ 4 . However, it is 
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possible to utilize a frequency lock-in amplifier to boost this limit to 10 -4 — 10~ 6 . For the sake 
of argument, we set a conservative number of 10 -4 as a practical detection limit. Considering 
a symmetric set-up such that o~\ = o-i = cr, from equation 3 one immediately finds that 021/cr 
is on the order of t^-ct, with additional factors determined by the ratio between the Green 
function and local density of states. For a case of ballistic transport through surface states 
(see below), \G\ 2 / p 2 is approximately f^, where k is the wave number at energy p, and r is 
the distance between the tips. The closest tip-to-tip distance is determined by the radius of 
curvature and the aspect ratio of the tips. Recent advances of tip-fabrication techniques can 
reproducibly make high aspect ratio tips with a radius of curvature on the order of 50 - 100 
rA [fF|. It is thus conceivable to consider operating a double-tip STM in the range of 300 - 
1000 rA tip-to-tip distance. This gives \G\ 2 /p 2 on the order of one percent. Since h/e 2 = 25 
Ml, one can immediately identify the practical operation range of tunnel junction resistance 
to be on the order of 1-10 MO. For the case of an anisotropic Fermi surface (as discussed 
below), \G\ 2 /p 2 is on the order of unity in the same range of consideration and the STM 
junction resistance can be as large as 1 GQ. The most advantageous case would be a one 
dimensional system, such as the fullerene nanotube, where the transconductance does not 
depend on distance, and is of order unity. In the following, we describe some of the most 
important applications of the double-tip STM. 

-Ballistic Transport and Surface State Band Structure. Consider a situation where 
surface states play a dominant role in electron transport. It is well known that the Green 
function for 2D free electrons is given by 

im TT m n , im I 2 i(kr _E) , A , 

2h 2 V ; 2h 2 V nkr y ' 

where k is the wave number, r = |r2 — T\\ is the tip-tip distance, and Hq is a Hankel function 
of the first kind. The right hand side is valid in the asymptotic region kr » 1, which is also 
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the regime of our interest. Therefore, the transconductance decreases inversely with r and is 
isotropic. On an actual crystalline surface, the states are Bloch waves e jkr u n fc(r) with band 
energies e n (k). It can be shown that for large kr, the Green function can be approximated 
as [§ 

G(r a , r 2 ; e) « -m nkc ( ri )< kc (r 2 ) e^^-O-f] (5) 

where ku and fc^ are the components of k in the directions parallel and perpendicular to 
r 2 — ri, respectively. The partial derivatives are evaluated at the point k c on the energy 
surface where the normal vector (group velocity) points in the direction of r 2 — ri. In 
the following analysis we assume for simplicity that there is only a single or single dominant 

Like the free electron case, the transconductance has an overall inverse r dependence. 
However, two interesting new features appear for a crystal: (1) The transconductance is 
modulated by the magnitude squared of the Bloch functions. (2) There is also an overall 
orientational dependence from the factors involving the partial derivatives of the band en- 
ergy. In Fig.2, the transconductance for a square lattice is plotted at an energy close to a 
nested energy surface, showing an inverse r dependence and a pronounced anisotropy. Bloch 
function modulation has been removed by averaging over a unit cell. Since the critical point 
k c runs through the energy surface as one changes the orientation of r 2 — r 1; it is possible 
to reconstruct the constant energy contour of the surface band structure for the filled and 
empty states using the mapped-out |G| 2 . The oscillatory modulation of the transconduc- 
tance should also tell us the shape of the Bloch waves for each k c . The energy resolution 
of the tunneling measurement is practically limited only by ksT. Sub-meV resolution is 
routinely obtained ||. In comparison, angle- resolved photoemission (ARPES) can map-out 
k-dependent band structure only for the filled states and its resolution is currently limited to 
about 15 meV both by the spectrometer and the photon source. A fourth power dependence 
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of the signal-to-noise ratio on the demanded resolution makes it very difficult for ARPES to 
achieve a finer energy resolution. 

-Phase shifts from a surface defect. One can move the tips near (but still far compared 
to the wave length) to a surface defect and observe how it scatters the electrons by an 
interference effect. An electron may propagate freely from tipl to tip2, or it may propagate 
to the defect and be scattered to tip2. The superposition of these processes can give rise to a 
modulation of |G| 2 , when tip2 is moved around relative to tipl and the defect. To illustrate 
the point, consider the 2D free electron model again. If the scattering is dominated by 
the s-wave channel, the interference pattern will consist of curves of constant path-length- 
difference 7"i + r 2 — r, where is the distance from the defect to tip-i. Clearly, these are 
hyperbolic curves, with the positions of tipl and the defects being the two focal points. 
Quantitatively, we have 

\G\ 2 /\G Q \ 2 ^l + J-^—wi9ca&{k{r-T 1 -r 2 )-'^-9) (6) 
V 7r£rir 2 4 

where 9 is the phase shift of the defect. The phase shift can thus be determined by the 
strength and the position of the modulation. 

In the presence of additional scattering channels, the above formula is generalized as 

/ 8^ 00 7f 

|C| 2 /|C 1 2 ~ 1 + \ — Yl cos(ma) sin# m cos(Hr - r\ - r 2 ) - mn - — - 9 m ) (7) 

V *ir 2 4 

where 9 m is the phase shift for the mth angular momentum. Now, the interference pattern 
also depends explicitly on the angle a between the vectors from the defect center to the tips. 

On an actual solid surface, where defects scatter Bloch waves, the above results should 
remain valid as long as the energy surface is fairly isotropic. In the general situation, one 
has to use the scattering matrix Tk cl k c2 , where k cl is a Bloch wave vector corresponding to 
a group velocity in the direction from tipl to the defect, and k c2 corresponds to a group 
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velocity in the direction from the defect to tip2. The transconductance should behave like 

\G\ 2 = |Go(r 1 ,r 2 ) + G (r 1 ,0)T kcl!kc2 G (0,r 2 )| 2 (8) 

where Go is the Green function in the absence of scattering. The asymptotic form of this 
result at large distances can be obtained by using (5). 

A single-tip STM can only see defects having a substantial disturbance of the local 
density of states on a sample surface. The double-tip STM should be able to detect the 
'mines' buried fairly deep under a surface though the 'radar' of scattering interference. Take 
again the free electron model (now 3D) to illustrate the point. In terms of the positions of 
tip2, the contours of constant path length difference are now the intersections of the sample 
surface with hyperbolic surfaces whose focal points are at tipl and the defect center. The 
defect can then be located by the characteristics of the interference pattern, e.g. it must lie 
under the symmetry axis of the pattern. 

- Transition from Ballistic to diffusive transport to localization. In the absence of defects, 
the transconductance measures ballistic electron transport between the tips. As seen above, 
for surface states, this is signified by an inverse r dependence, modulated by an angular 
dependence if the energy surface is not circular. When a defect is present in the neighborhood 
of the two tips, scattering interference will occur. The interference pattern will become more 
and more complicated if more defects are included in the way between the tips. When the 
tips are moved apart by a few elastic mean free paths, diffusive electron transport should 
begin to take place. It can be shown that |G(ri, r 2 ; e)| 2 measures the time integral of the 



probability for a wave packet of energy e to move from to r 2 in a given time t ||10|| . In the 
diffusive regime this probability goes as 

2 

(7r J Dt)- 1 exp(-^) (9) 

where D is the diffusion constant. Here, we should emphasize that this formula only describes 
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the overall trend; statistical fluctuations can still persist in the diffusive regime in the form 



of the so called universal conductance fluctuations [11]. The average transconductance in 
this regime should then behave as 

(nD^E^/DQ (10) 

where Ei(X) = ^-e~ x is the incomplete T function, and t c is a cutoff time beyond which 
elastic diffusive behavior ceases to occur. In the absence of inelastic scattering, this cutoff 



time is given by l^/D, where l c is the Anderson localization length JTT[. It is interesting to see 
that the transconductance for r << l c behaves like (irD)~ 1 2 ln(/ c /r). This slow fall off with 
distance is in sharp contrast with the behavior in the ballistic regime. Also, the angular 
dependence should go away for a surface with square or hexagonal crystalline symmetry, 
because, as a second rank tensor, the diffusion coefficient cannot distinguish such point 
group symmetries from full isotropy. Finally, when r is beyond the localization length, the 
Green function and the transconductance should drop exponentially. 

-Inelastic mean free path. In the above discussions, we have ignored electron-electron 
and electron-phonon interactions in order to simplify the presentation. It is well known 
that the single particle Green function will acquire a self energy with an imaginary part as 



well as a real part due to such interactions | 12fl . The imaginary part pushes the poles of 



the Green function off the real axis, yielding an exponential decay of the Green function 
in distance r = |r 2 — r i|- Physically, an electron tunneled in from tipl may lose energy 
by exciting electron-hole pairs or emitting phonons as it travels in the sample, and become 
unable to tunnel out to tip2. The typical length scale over which such a process occurs 
defines the inelastic mean free path, and is given by the decay length of the Green function. 
Therefore, the inelastic mean free path and its energy dependence may be measured by 
observing how the transconductance decays with tip-tip distance, and how the decay length 
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varies with \i — /i . For /i — /i above the Debye energy (tens of meV, hot electrons), phonon 
emission dominates the inelastic processes. At lower energy differences, electron-hole pair 
excitations become dominant. In a Fermi liquid without disorder, the mean free path goes 
as hE F vp/(n — /io) 2 - A different energy dependence has been predicted for non-Fermi liquid 
systems. The ability to measure the energy-dependent mean free path using the double-tip 
method should have important impact on this issue 

-Summary and Discussion. Because a double-tip STM can probe the all important 
single particle Green function of a sample, it has the potential of becoming an extremely 
useful new tool in surface analysis. We have identified key experimental parameters for such 
measurements. We have further described some basic applications of a double-tip STM: (1) 
Probing k resolved band structure of surface states, as well as the shape of Bloch functions; 
(2) Measuring scattering phase shifts or amplitude of surface defects; (3) Observing transition 
from ballistic to diffusive transport to localization, and (4) Measuring inelastic mean free 
paths. 

Apart from these applications, one can consider applying a magnetic field to observe the 



cyclotronic motion of electrons in the semiclassical regime [13|, or to possibly probe some 



properties of a quantum Hall system in a strong magnetic field |]14| . One can also consider 
mapping out the quasi-particle band structure of a superconductor through ballistic transport 
or gap structure through Andreev reflections as proposed in Ref.0. 
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Figure Captions: 

Fig.l Schematic diagram of the double-tip STM experimental set-up. Tip 1 is biased at 
V\ and tip 2 at V 2 relative to the sample, h and I 2 are also measured relative to the sample. 
When Hi > fi 2 , h contains a transconductance component due to the co-tunneling process. 

Fig. 2 The distance and angular dependence of the transconductance for a surface with 
square symmetry. X and Y are the coordinates of r 2 — ri (100 ~ 600 lattice constants) on 
a sample surface. The vertical axis is the square of the Green function normalized by the 
local density of states, which is the same as the transconductance normalized by the junction 
conductances up to a constant factor. The corresponding energy surface (nearly nested) is 
shown on top of the figure. 
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Tip 1 at Vj 



Tip 2 at V 2 




